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Abstract 

Existing algorithms for isolating real solutions of zero-dimensional polynomial 
systems do not compute the multiplicities of the solutions. In this paper, we de- 
fine in a natural way the multiplicity of solutions of zero-dimensional triangular 
polynomial systems and prove that our definition is equivalent to the classical def- 
inition of local (intersection) multiplicity. Then we present an effective and com- 
plete algorithm for isolating real solutions with multiplicities of zero-dimensional 
triangular polynomial systems using our definition. The algorithm is based on in- 
terval arithmetic and square-free factorization of polynomials with real algebraic 
coefficients. The computational results on some examples from the literature are 
presented. 

1 Introduction 

Real solution isolation for polynomials/zero-dimensional polynomial systems/semi- 
algebraic systems is one of the central topics in computational real algebra and com- 
putational real algebraic geometry, which has many applications in various problems 
with different backgrounds. 

The so-called real root/zero/solution isolation of a polynomial/zero-dimensional 
polynomial system/semi-algebraic system with k distinct real solutions is to compute 
k disjoint intervals/"boxes" containing the k solutions, respectively. To our knowledge, 
designing algorithms for real root isolation for polynomials with rational coefficients 
was initiated by |4| in 1976, which was closely related to the implementation of CAD 
algorithm |3|. Designing and implementation of such algorithms have been deeply 
developed by many subsequent work ||5l [T] \W\ [TT] [TJl since then. Those algorithms 
are mainly based on Descartes' rule of sign or Vincent's theorem. 

To generalize the algorithms for polynomials to zero-dimensional triangular poly- 
nomial systems, one must consider real root isolation for polynomials with real alge- 
braic coefficients. There are indeed some work to generalize Descartes' rule of sign to 
polynomials with algebraic coefficients. However, dealing with algebraic coefficients 
directly may affect efficiency greatly. 

In ifm ITSl we considered real solution isolation for semi-algebraic systems with 
finite solutions. We introduced a method which always enables us to avoid handling 
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directly polynomials with algebraic coefficients and to deal with polynomials with ra- 
tional coefficients only. A recent algorithm in (T\ can compute the parity of the solu- 
tions as well as isolate real roots of zero-dimensional triangular polynomial systems. 

In this paper, we define in a natural way the multiplicity of solutions of zero- 
dimensional triangular polynomial systems and prove that our definition is equiva- 
lent to the classical definition of local (intersection) multiplicity. Then we present 
an effective and complete algorithm for isolating real solutions with multiplicities of 
zero-dimensional triangular polynomial systems using our definition. The algorithm is 
based on square-free factorization of polynomials with real algebraic coefficients and 
our previous work ITSl . We also provide computational results on some examples from 
the literature. 

In this paper, all polynomials are in C[X] = C[a;i, . . . , Xn] if not specified. 



First, let's recall the definition of local (intersection) multiplicity. We follow the nota- 
tions in Chapter 4 of |6|. Although some notations and definitions can be stated in a 
more general way, we restrict ourselves to the ring C[X] = C[xi, . . . , Xn] since we are 
interested in the complex or real zeros of zero-dimensional polynomial systems. 

For p = {t]i, . . . , rjn) € C", we denote by Mp the maximal ideal generated by 
{xi — rji, . . . ,Xn — 77,1} in C[X], and write 



It is well-known that C[X]mp is the so-called local ring. 

Definition 1 ||6| Suppose I is a zero-dimensional ideal in C[X] and p G Zero(/), the 
zero set of I in C. Then the multiplicity ofp as a point in Zero(/) is defined to be 



That is, the multiplicity ofp is the dimension of the quotient space <C[X]mj, / I'C[X]mj, 
as a vector space over C. 

For a zero of a zero-dimensional triangular set, there can be a natural and intuitional 
definition of multiplicity as follows. 

Definition 2 For a zero-dimensional triangular system. 



and one of its zeros, ^ = • ■ • , "^n), the multiplicity of S, is defined to be Y\ rrii, 

i=l 

where rrii is the multiplicity of Xi — ^ias a zero of the univariate polynomial fi{^i, . . . , Xi) 
fori^l,...,n. 



2 Multiplicities of zeros of triangular sets 




dimkC[X]MjIC[X]M,. 




n 
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Example 1 Consider the following triangular system: 



51 = + 2x1 + = 0' 

92= X2+x\+ XiX2 = 0, 
33 = x1+ XiX^ + XiX2 = 0. 



Let's compute the local multiplicity of (0,0,0) by Definition^ The multiplicity of 
Xi = 0, a zero ofgi, is 3. Substitute xi = in /2 and the resulted 52 w 92 = x^- Thus, 
the multiplicity of X2 = Q, a zero of g'2, is 2. Finally, substitute xi — X2 = in 53, 
and the resulted 53 is g'^ — x§. Thus, the multiplicity of x^ ~ 0, a zero of g'^^ is 2. As a 
result, the local multiplicity of (0, 0, 0) w 3 x 2 x 2 = 12. 

In the following, we will prove that Definition|2]is equivalent to Definition[T] Many 
notations and results are taken from ||6l. 

Usually, a total order that is compatible with multiplication and that satisfies I > Xi 
for all i's, is called a local order. 

Definition 3 / 9 /(Negative lexicographical ordering j 

Assume a = {ai, . . . , a„) e and (3 ^ {(3i, . . . , /3„) e We say X" >„; ;/ 
3i (1 < i < n) A (Vj, 1 < j < i =^ aj = (3j) A {a^ < Pi). 

Remark 1 The negative lexicographical ordering >„; is obviously a local order 

For a given order, lc(/),lm(/) and lt(/) denote the leading coefficient, leading 
monomial and leading term of /, respectively. For a set 5, It(S') = {lt(/) : / G 5}. 

Definition 4 [6| Lef R — C[X]m ant/ 1 (Z R be an ideal. A set {51, . . . , g„i} C / 
called a standard basis /or / with respect to <„; ;/ (lt(/)} = (lt((7i), . . . , \t{gm)). 

For a — (ai, . . . , «„) G N", define \a\ = X^i*^*- Po'" ™y polynomial g = 
CaX"" e C[X] with total degree d, we will write g'' = J2a Ca^''"'"'^" for the 
homogenization of g with respect to t. 



Definition 5 M Define f'X"' >'^i t^X^ if a + |a| > 6 + \[3\ ora^\a\=b+ but 

X" >nl X^ 



It is easy to verify that >',j; is a monomial order over C[i, X\. 
Theorem 1 ^^{Analog ofBuchberger's Criterion) 

Let G = {51, . . . ,grn}, > be any local order, and I be the ideal in C[X]mp gen- 
erated by G. G is a standard basis for I if and only if applying Mora normal form 
algorithm to each S-polynomial formed from elements of the set of homogenizations 
G'^ ~ ^g\, . . . ^ 5^} yields a zero remainder. 

For our purpose, we state the criterion in another form as follows. 

Theorem 2 Let notations be as in Theorem [7] G is a standard basis if and only if 
for any nonzero S-polynomial of g^ and g'j, denoted by Sij, there exist homogeneous 
polynomials U, Ai, . . . , Am G C[t, X] such that 




(1) 
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where lt(J7) = t"^ for some a, 

a + deg(5y ) = deg{Ai) + deg{g\') 
for all I whenever Ai ^ 0, and\t{Aig'^) <[^i lt(C/5'y ). 

Remark 2 We omit the proof of Theorem^ which is almost the same as that ofTheo- 
remUl The criterion in Theorem\2\is independent to any algorithms. One can use Mora 
normal form algorithm to get such representation as (O/or each Sij ifG is a standard 
basis. 

Without loss of generaUty, in the rest of this section we assume p = (0, . . . , 0) is a 
zero of the triangular set under discussion and focus on its multiplicity. Consider the 
following triangular set with leading terms cix'"^ , • • ■ , Cn^:^™" respectively w.r.t. the 
order >„;: 

{fi{xi) = cix'"' + ti{xi), 
f2ixi,X2) = C2X^^ +t2ixi,X2), 
fn {-^l : • ■ • : -^71 ) — ^n-^n ^ ~^ ('^l i ■ ■ • 5 -^n ) ; 

where ti{xi, . . . , Xi) is a polynomial in xi, . . . , for i — 1, . . . , n and q's are con- 
stants. Without loss of generality, we assume q's are all 1 in the proof of the following 
proposition. 

Proposition 1 Let T be as above and I — (T) the ideal generated by T in the local 
ring C[X]^2:i....,x„)- Then T is a standard basis for I with respect to >„;. 

Proof. According to Theorem |2] we only need to show that every nonzero S- 
polynomial of each pair of T'^ = j/^*, . . . , /^} can be represented in the form of 

Assume that = + /j' = t''x"^' + tj and a < b. The S-polynomial, 

of//' and /j' is 

C _ 4-b-a inj rh _ rrii rh 
tJi] — i' -^j Ji -^i Jj ■ 

Let pi — t^^°-x^' fl^ and p2 — a;™'//- Under the order <'^;, the first term of pi is 
equal to the first term of p2. If Sij ^ 0, there exists some L such that under the order 
<'^i the L-th term of pi is not equal to the L-th term of p2 and the fc-th term of pi is 
equal to the /c-th term of p2 for all 1 < fc < L. Then, the fc-th terms of pi and p2 can 
be represented as x"^^qk and t^~°'x^^qk for some qt, respectively. Thus, and fj- 
can be respectively rewritten as 

= xT' (i" + Q) + T2, = t^-'-x^' {t" + Q)+7i2. 

L-l 

where Q = ^ q^ and /j2 and fj2 are the remained parts of and respectively, 

fc=i 

which satisfy that 

Ht'-^x';'jr2)^HxrT2)- 

It is easy to verify that 5^ = t^^^x'^'J^ - xf^J~^. Then 

lm(5y) = max(lm(i''-"x"^i^;),lm(a:™'7^)) 
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under the order <^j. Thus, 

Let U ^t" + Q, Aj = Tv2,Ai^J~^ and Ai =Q{l^i and I ^ j). Then we have 

n 

1=1 

and all the requirements in Theorem |2] are met. Thus T is a standard basis for {T) w.r.t 

<nl- 

In order to prove the equivalence of Definitions |2] and [T] about the (local) multiplic- 
ity, we need the following theorem, which can be found in ||6l. 

Theorem 3 Let I be an ideal in a local ring R, and assume that dim^ R/ (lt(/)} 
is finite for some local order >. Then we have 

dimfc i?// = dimfe i?/(lt(/)). 
Theorem 4 Let notations be as above and T a zero-dimensional triangular set with a 

n 

zero p = (0, . . . , 0). If the multiplicity of p defined by Definition^is rn = Y[ then 

i=l 

the local multiplicity, defined by Definition\l\ of p as a point of (T) is also m. 

n 

Proof. If the multiplicity is Y\ "^j the sense of Definition |2l then T can be 

i=i 

rewritten as 



T 



{ /i(xi) = (ci+tn(xi)X\ 

f2{Xl,X2) = (C2 +t2l{xi,X2))x^^ + Xit22{xi , X2) , 
. . . , 

n-1 

fn(X) = (C„ +<„i(^)X'" + E X.^r^r+l{X), 

i=l 



where X = (xi, . . . , Xn), the (X)s are polynomials in (xi, . . . , Xi) and the tii{X)s 
do not contain constants. 

Under the order >„;, the leading term of fi{xi, . . . ,Xi) is Cix"^^ for i — 1, . . . ,n. 
According to Proposition[Tl T is a standard basis of / = (T) . Thus, 

(lt(/)) = «\...,C">- 
Let R = CfXj^j-j 3;^^. According to Theorem[3] 



dimfe R/IR ^ diuik 1 , . . . , x^^" )^ = H 



3 Algorithm for real solution isolation with multiplicity 

In this section, based on the results in last section, we present an algorithm for real root 
isolation with multiplicity of zero-dimensional triangular polynomial equations. That 
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is, we not only isolate the real roots, but also compute the multiplicity of each real root 
by Definition |2] at the same time. In this section, the input polynomial or polynomial 
set to our algorithms is taken from Q[X]. 

It is well-known that there exist some efficient algorithms for real root isolation of 
polynomials or polynomial equations or semi-algebraic systems P1l5l[n [T0lll7|[m i2l 
[181 . To obtain the multiplicities of the real roots at the same time, our idea is simple 
that is to take use of square-free factorization of polynomials with rational or algebraic 
coefficients. When dealing with algebraic coefficients, we make use of the idea in 
ifTTlfTSl which enables us to deal with rational coefficients instead. 

For the univariate case, suppose p — HiLi Pi- Isolating the real zeros of p with 
multiplicity contains two main steps. One is to compute the squarefree factorization of 
p, the other is to isolate the real zeros of the squarefree part of p. We can use many 
existing tools to obtain the squarefree factorization, i.e., those p^s. Then we know at 
once the multiplicities of those real zeros of each pi. In principle, we may isolate the 
real zeros of the squarefree part of p in two ways. One way is to isolate the real zeros 
of P1P2 ■ ■ - Pk first and then match the zeros with pi to obtain correct multiplicities. 
The other way is to isolate the real zeros of each pi separately. However, in the later 
way, we may need to compute a root gap of p first. Anyway, the univariate case can be 
efficiently dealt with. So, we do not enter the details of such algorithms and only give 
a description of the input and output of such function. 

Calling sequence UniIsol(/(x)) 
Input: a univariate polynomial /(x) 

Output: a set of elements of the form ([a, b],m) where [a, b] is an interval containing 
exact one real root of f{x) = and m is the multiplicity of the root. There are not any 
real roots of f{x) — outside the intervals. 

Then let us consider the multivariate case. To be more precise, we state our problem 
as follows: that is to isolate the real solutions with multiplicities of the following zero- 
dimensional triangular polynomial set 

T = {/l(a;i), f2{xi,X2), fn{xi, . . . , Xn)}. 

In principle. Definition |2] suggests a naive method to compute the local multiplicity as 
follows. First compute all the zeros of fi{xi) and their multiplicities by Unilsol; then 
"substitute" the zeros for xi in f2{xi, X2) one by one, and compute all the zeros of the 
resulted f2{xi, X2) and their multiplicities by Unilsol again, and so on. Of course, 
in general we cannot directly substitute the zeros in those polynomials because they 
may be algebraic numbers of high degrees. Nevertheless, this naive method is the main 
framework of our algorithm. 

LetT, = {fi{xi)j2{xi,X2),...,fi{xi,...,Xi)}. We will call 

{[ai,bi],...,[ai,bi]) or ([ai, 5i], . . . , [a^, 6^], to) 

an interval solution of Ti (with multiplicity m) if the "box" [ai, 61] x • ■ • x [a.;, bi] 
contains exact one real solution of Ti (and to is the multiplicity of the solution). 
If Ti has k distinct real solutions, a set of k interval solutions of Ti containing re- 
spectively the k real solutions is called a solution set of Ti. For an interval solution 

r = ([ai,5i], . . . , [ai,bi]), we define Nr = [xi ~ ai,6i ~ xx, . . . ,Xi - ai,bi - x^ 
and A^r > stands for a 1 < x\ < &i,...,a,; < Xi < bi, i.e., {xi,...,Xi) £ 
[ai,6i] X • • • X [ai,bi]. 
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Suppose we already have a solution set of Ti and 

(6, • ■ • ,6) e [ai,bi] X • • • X [ai,bi] 

is a real root of Ti with multiplicity m. To isolate the real zeros of ■ ■ ■ Xi+i) 

with multiplicity, we need to 

1. compute the algebraic squarefree factorization of ■ ■ ■ Xi+i), and 

2. isolate the real zeros of the squarefree part computed. 

Let us first consider the second task, i.e., how to isolate the real zeros of /i+i i , . . . , 
^i, Xi+i) if it is squarefree. In 1 18 1, we proposed a complete algorithm, called RealZeros, 
for isolating the real solutions (without multiplicities) of semi-algebraic systems. Our 
second task can be accomplished by a sub-algorithm of RealZeros. The key idea 
of the algorithm is to compute two suitable polynomials fi+i and /i+i in Xi+i with 
rational coefficients such that 



fi+l < fi+l{^l, • ■ • I < fi+i 

by using interval arithmetic and those intervals [ai, . . . , [a^, bi\. And the real zeros 
of . . . , ^i, Xi+i) can be isolated through isolating the real zeros of /i+i and 

fi+i- Therefore, we can avoid dealing with polynomials with algebraic coefficients 
directly. In the following, we call this sub-algorithm Algebraiclsolate. 

Calling sequence AlgebraicIsolate((ji(a;i, . . . ,Xi+i),Ti,r) 
Input: a squarefree polynomial g{xi, . . . , Xi+i), a zero-dimensional triangular poly- 
nomial set Ti as above and an interval solution r = {[ai, bi], . . . , [ui, bi]) which con- 
tains exact one real zero {£,1 , ■ ■ ■ , of T^. 

Output: a list of isolating intervals of real zeros of 5(^1, . . . , ^i, x^+i). 

For the detail of the algorithm Algebraiclsolate, please be referred to ifTsl . 

Now, we turn to the first task, i.e., compute the algebraic squarefree factorization 
of . . . ,£i,Xi+i). One may use some existing algorithms for algebraic factor- 

ization, see for example llT4l . to accomplish the task. In the following, we propose a 
method for algebraic squarefree factorization based on algebraic gcd computation. A 
key manipulation in the computation is to count real solutions of semi-algebraic sys- 
tems by an algorithm RealrootCountingin lfT6lFl 

Calling sequence RealrootCounting(_F, N, P, H) 

Input: a zero-dimensional polynomial set F, a list of non-strict inequalities N , a list 
of strict inequalities P and a list of inequations H. 

Output: the number of real roots of the system {F = 0, > 0, P > 0, ^ 0}. 

Calling sequence AlgebraicGCD(pi(a;i, . . . ,Xi+i),p2{xi, . . . ,Xi+i),Ti,r) 
Input: two polynomials pi , p2 in xi , . . . , Xi+i, a zero-dimensional triangular polyno- 
mial set Ti as above and an interval solution r — ([oi, 61], ... , [ai, bi]) of Ti. 
Output: The greatest common divisor of pi and p2 viewed as polynomials in Xi+i 
w.r.t. the interval solution r, i.e., gcd(pi(^i, . . . , Ci, Xi+i),P2{£.i, £,i,Xi+i)) where 
(^1, . . . , ^i) is the only real solution in r. 

'The algorithm is called nearsolve in ^16] 



7 



Step Suppose the subresultant chain of pi and p2 w.r.t. .t,;+i is S'p, . . . , Sq 

with principal subresultant coefficients . . . , i?o, respectively. Setj ^ 

0. 

Step 1 Compute Rj. 

Step 2 If RealrootCounting(ri, Nr, [ ], [Rj]) = 0, i.e., the interval solution makes 
Rj vanish, then set j ^ j + 1 and go to Step 1 . 

Step 3 Return 5*^ . 

There are several mature algorithms IS] [T3l for squarefree factorization of polyno- 
mials in K[x] where IK is Z, Q or a finite field. It is well known that such algorithms 
for univariate case only contain two main manipulation: gcd computation and poly- 
nomial division in K[x]. If we replace the gcd computation in those algorithms with 
our AlgebraicGCD computation and replace the division manipulation with pseudo- 
division, then those algorithms will compute algebraic squarefree factorization as we 
want. So, We only give a simple description of our algorithm here. 

Calling sequence AlgebraicSqf reeFactor(p(2;i, . . . , Xi+i), T^, r) 
Input: a polynomial p in xi, . . . , Xi+i, a zero-dimensional triangular polynomial set 
Ti as above and an interval solution r — {[ai, bi], . . . , [ui, bi]) of T^. 
Output: the squarefree factorization of p viewed as a polynomial in Xi+i w.rt. the in- 
terval solution r, i.e., the squarefree factorization of p{^i , . . . , , Xi+i) where {^i, ■ ■ ■ ,£,i) 
is the only real solution in r. 

Now, we are ready to describe our algorithm Multilsolate for real solution iso- 
lation with multiplicity of zero-dimensional triangular polynomial sets. 

Calling sequence Multilsolate(T) 

Input: a zero-dimensional triangular polynomial set T = {/i(a;i), . . . , fn{xi, ■ ■ ■ ,a;„)}. 
Output: a solution set of T with multiplicity. 

Step 1 i^l,Li^ UniIsol(/i). 

Step 2 Li is a solution set of T,; with multiplicity. If i — n, return i„. 

Step 3 For each interval solution r = ([oi, . . . , [a^, bi]) in Li with multiplicity, 
compute AlgebraicSqf reeFactor(/i_|_i(xi, . . . , ^i+i), Ti, r). So, we know 
at once the multiplicity of each factor Assume is the squarefree part of 
fi+i- Then, by applying AlgebraicIsolate(/,;+i, T^, r) we can obtain the 
isolating intervals of real zeros of fi+i- So, it is easy to obtain a solution set 
Li+i of Ti+i with multiplicity by Definition|2] 

i ^ i + 1 and go to Step 2. 

Remark 3 Let r ~ ([ai, . . . , [a„, 6„]) be an interval solution of T and ^ = 
(Cij ■ • ■ iCn) the real solution in r. ^lc(/i)(^i, . . . 7^ 0/or 2 < i < n, T 

is said to be regular w.r.t. ^ (or r). If fiiXi, . . . , ■Ci-ii ^i) is squarefree for I < i < n, 
T is said to be squarefree w.rt. ^(orr). 

It is clear that Multilsolate(T) actually computes as well a regular and square- 
free decomposition of the given triangular set T w.r.t. its real zeros, respectively. That 
is to say, we compute a set of triangular sets Wj and their solution sets Qj such that 
UjQj is a solution set ofT and each Wj is regular and squarefree w.rt. each solution 
in Qj. If we modify slightly the algorithm, we can output the regular and squarefree 
decomposition. 
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4 Examples 



The algorithm Multilsolate has been implemented as a Maple program which is 
included in our package DISCOVERER {15\. For an input zero-dimensional triangu- 
lar system, our program can compute the real solution isolation of the system with 
multiplicity and output a regular and squarefree decomposition (see Remark[3]l of the 
system w.r.t. those real solutions. Our program can detect whether the input system is 
zero-dimensional. If it is not, the program will return a message: "The dimension of 
the system is positive." 

In this section, we illustrate the function of our program by some examples. The 
timings are collected on a Thinkpad X200 running Maple 11 with IG memory and 
Windows Vista. 

Example 2 Consider the following triangular system: 

/i = 2; - 2, 

f2 = ix + y-3f{y + 3), 
fs = {yz^ + XZ + - y^z + x~y). 

Within 1.6 seconds, our program outputs a solution set as follows. 

[ [ [[2, 2], [-3, -3], [-i 0]], 1 ], [ [[2, 2], [-3, -3], [1, 1]], 2 ], 

[ [[2, 2], ha, -3], [-1 -i]], 2 ], [ [[2, 2], [1, 1], [-1, -1]], 15 ] ]. 

That means the system has 4 real solutions which are of multiplicities 1, 2^2, Ib^ re- 
spectively. Our program also outputs a regular and squarefree decomposition of the 
system w.r.t. the four distinct real solutions respectively as follows. 

[x-2,?/ + 3,l + 125z], 
[a;-2,?/ + 3,-l + 3z2-2z], 
[x~2,y^l,z + l]. 

Note that the second and third solutions are both solutions to the second equations 
above. 

Example 3 Consider the following triangular system: 

/i = (x + l)(x-2), 
/2 = {x-y + lf{y-^) + {y-3)x, 
/3 = (a;j/-6)z2 + 2z + l. 

The system has 7 real solutions all of multiplicities 1. The computation costs 0.7 sec- 
onds. 

[ [ [[2, 2], [3, 3], [-1/2, -1/2]], 1 ], [ [[-1, -1], [-1, -3/4], [-3/8, -1/8]], 1 ], 
[ [[-1, -1], [-1, -3/4], [3/8, 7/8]], 1 ], [ [[-1, -1], [1/2, 3/4], [-3/8, -1/8]], 1 ]. 
[ [[-1, -1], [1/2, 3/4], [3/8, 3/4]], 1 ], [ [[-1, -1], [5, 21/4], [-3/8, -1/8]], 1 ], 
[ [[-1,-1], [5, 21/4], [1/4, 1/2]], 1] ]. 

A regular and squarefree decomposition is 

[a;-2,/,l + 2z], [x+lj,g\, 
where f = x^y—5x^—2xy'^ + 13xy—13x+y^ — 7y'^ + lly — 5, g — yz^x—6z'^+2z+l. 
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Example 4 The following triangular system is taken from ^7^. 



Within 0.1 seconds, we obtain two distinct real roots with multiplicities 16. 



And a regular and squarefree decomposition is 



Example 5 The following triangular system is taken from 



\ /2 = 2/^ + xy^ + 3j/2 - 6x^2/2 + 4a;?/ + 2xy'^ - Ax'^y + 4x + 2. 
The time for computation is 3.6 seconds and we obtain 12 distinct real roots. 



[ [[51/32, 13/8], [-119/32, -475/128]],! ], 

[ [ [51/32, 13/8], [-147/128, -145/128] ],1 ], 
[ [ [51/32, 13/8], [53/64, 107/128] ],1 ], 
[ [[51/32, 13/8], [307/128, 77/32] ],1 ], 
[ [[-5/8, -19/32], [-3/8, 1/4] ],1 ], 
[ [[-5/8, -19/32], [13/8, 17/8] ],1 ], 
[ [ [45/32, 23/16], [-3025/1024, -1499/512] ], 1 
[ [ [45/32, 23/16], [-1347/1024, -2639/2048] ], 1 
[ [[45/32, 23/16], [11/8, 3/2]], 2 ], 
[ [[-23/16, -45/32], [-5/8, -1/8] ],1 ], 

[ [[-23/16, -45/32], [17/4, 5] ],1 ], 
[ [[-23/16, -45/32], [-3/2, -11/8]], 2 ]. 



It is clear that two of the solutions are of multiplicities 2 and the others are of multiplic- 
ities 1. With respect to those solutions, we have a regular and squarefree decomposition 
as follows. 



[x^-x-\,hi\, [x^-2,h2], [x2-2,ft3], 



hi^y^ + xy^ + 3y'^ - 6x^2/^ + Axy + 2xy'^ - Ax^y + 4x + 2, 

h2 = -23354573041809- 9122537689096x7/2 + 39406733143725xy + 

17148617740054X+ 13135577714575y2 - 54735226134576j/, 
/i3 = -104xy + 335y - 335x + 208. 




[ [[0,0], [0,0], [-1,-1]], 16 



[ [[0,0], [0,0], [0,0]], 16 ]. 



[x,y, z + - 7x^ - 8x2]. 




where 
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